Cosmic String Lens Phenomenology: Model of Poisson Energy Distribution 
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We present a novel approach for investigating lens phenomenology of cosmic strings in order to 
elaborate detection strategies in galaxy deep field images. To account for the complexity of the 
projected energy distribution of string networks we assume their lens effects to be similar to those 
of a straight string carrying a random lineic energy distribution. In such a model we show that, 
unlike the case of uniform strings, critical phenomena naturally appear. We explore the properties 
of the critical lines and caustics. In particular, assuming that the energy coherence length along 
the string is much smaller than the observation scale, we succeeded in computing the total length 
of critical lines per unit string length and found it to be 4/\/3E(3/4). The length of the associated 
caustic lines can also be computed to be 16/(7r \/3) E(3/4). The picture we obtain here for the 
phenomenology of cosmic string detection is clearly at variance with common lore. 
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I. INTRODUCTION 

Although the current results on large-scale structure 
formation favor models with initial aiabatic scalar per- 
turbation [Q, the formation of cosmic strings is general 
enough to merit further investigations j^]. It is in par- 
ticular not excluded that strings may form after an in- 
flationary phase We are interested here in the phe- 
nomenological aspects of cosmic strings for lens distortion 
effects whereas most of the previous investigations have 
been focused on string detection from multiple quasar im- 
ages or Cosmic Microwave Background fluctuations 
With the advance of new generation of large CCD 
cameras the best direct evidence for cosmic string relics 
is however likely to be obtained from distortion effects 
they induce on background objects such as galaxies. 

In a companion paper we insist on the generic prop- 
erties expected for cosmic strings as far as lens effects 
are concerned: string lens effects are equivalent to those 
induced by a lineic energy distribution. This effective 
energy distribution should obviously take into account 
the lens energy density, its tension, as well as kinetic en- 
ergy that might result from rapid movements or energy 
currents along the string 0|. 

The explicit computation of string induced lens effects 
has been done in various cases, for long strings ||^, cir- 
cular loops in plane transverse to the line of sight [B|, 
or small loops with a multipole expansion approach 
However, much more complex situations have been ex- 
plored so far with numerical experiments only In 
particular the study presented in suggests a quite dif- 
ferent lens phenomenology for the images of background 
galaxies but lacks analytical insights. Having in mind 
such phenomenological effects, we try in this paper to 



get more quantitative results by introducing a simplified 
description of strings. In Sect. || we present the model 
of "Poisson" energy distribution we use. In Sect. Ill 
general phenomenological properties observed in such a 



model are presented. Sect. [V is devoted to more pre- 



cise calculations on the properties of the critical lines and 
caustics for such a model of strings and finally we evoque 
observational aspects. 



II. STRING MODEL 

A. General formalism 

In general lens effects are encoded in an amplification 
matrix. A, which describes the way the angular positions 
in the image plane are transformed to those in the source 
plane (see [|ll| for a comprehensive description of lens 
physics) . It is usually written in terms of the convergence 
field K and the shear field (71,72), 



A = 



1 



K - 71 

-72 1 



-72 
K + 71 



(1) 



that are directly related to the projected energy density. 
The remarkable property for lens effects are obtained 
when one of the eigenvalues of A is crossing or getting 
close to zero. The difficulty however is that these features 
are neither local nor linear in the energy density and lo- 
cation of the string. The fact that the projected density 
is likely to fluctuate substantially, because of local veloci- 
ties, wiggles, longitudinal motions etc., has somehow to 
be incorporated in the description. 



1 



In case of cosmic strings, the elements of the deforma- 
tion matrix can be derived the effective projected poten- 
tial (p{x, y) which can formally be written in terms of the 
angular positions (x, y) (e.g. Eq. (37) of as, 



(p(a;,y) = 4G 



D 



LS 



OS 



ds Ai[a;str.(s),2;str.(s)] X 



log {[x - Xst,Xs)f + [y~ yst,.{s)f) 



2\l/2 



(2) 



where (a;str.(s),2/str.(s)) are the angular string co- 
ordinates for the angular curvilinear position s, 
At(2^str. (s), 2/str. (s)) is the projected energy density at 
those positions, G is the Newton constant and Dls/-Cos 
is the ratio of the angular distance between the string 
and the source-plane to the one between the observer 
and the source plane in the thin leng|^ approximation 
I p^ . The projected energy density is a combination of 
the projected Tqo, Tqz and T^z components of the stress- 
energy tensor of the string if the line-of-sight is along the 
z direction. The displacement field is then given by 



£,i = -di(p{x,y) 



(3) 



and the elements of the deformation matrix can be writ- 
ten as. 



72 = 2dxdy (p{x,y), 



(4) 
(5) 



the local convergence being zero except on the string it- 
self. 



B. The Poisson string model 

It seems very difficult to take into account all the fea- 
tures that must be included: the string are far from be- 
ing straight lines with uniform energy distribution 
We choose to describe the energy fluctuation in a sim- 
ple manner, assuming that the string follows a straight 
line, but with local energy fluctuations. This fluctuations 
are assumed to account for the various changes of shape, 
density of the strings, for possible non-standard equation 
of states, or for the existence of currents along the string. 
We therefore assume the string to be straight along the 
y direction. 



a;str.(s)=0, ystr.(s) 



(6) 



and iJ,{s) to be a random field. Note that it does not mean 
however that the string is actually orthogonal to the line- 
of-sight. It simply means that the string (or rather the 



*The thin lens approximation is appropriate if one is inter- 
ested in only a fraction of a string spanning at most a few 
degrees on the sky. This would not be appropriate for appar- 
ent string crossings for instance. 



section of string we are interested in) appears roughly 
straight on the sky. To specify our model we still need 
to explicit the statistical properties of the fu, field. In 
general the results are going to depend on the chosen 
global statistical properties of /i and not only on its 2- 
point function for instance. In particular one question to 
ask is whether there is a finite correlation length along 
the string or not. Let us assume that the string can be 
described by the following properties. 



(m(s)) = Mo 
(m(si)m(s2))^ = C2Aio exp 
(^(si) . . ./x(sp))^ = CpM? X 
exp 



|Sl - S2| 



So 



max{sp} — min{sp} 



So 



(7) 

(8) 

(9) 



so that /i(s) has a finite coherence length, sq, and /i(s) is 
essentially uniform with typical value fj,Q at scales smaller 
than So- The statistical properties of the displacement 
field or of the elements of the amplification matrix are 
then described by the ones of /i. For instance, 



(7f(^)), 
dsi 
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+ 00 



dsp Gi(x, si) . . . Gi{x, Sp) X 



exp 
with 



max{sp} — min{sp} 



So 



Gi(a:,y)=4 



G2(x,y) = -8 



h x-' 
xy 



{y^ + x^) 



(10) 

(11) 
(12) 



The general expressions of these integrals are compli- 
cated, but they can be computed at leading order is sq/x. 
First of all we can notice that 



(71) = (72) = 0. 



(13) 



At leading order in sq/x the computation of the results 
can be obtained by an expansion of Gi{x, Sj) with respect 
to Sj in the vicinity of si (for instance) as 



Gi{x, Sj) 
It leads to 



(s, -si)f dPG, 



p 



dsP 



(x,si). 



(14) 



(7,^)«8c2G2 ^ls„2 
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dsiGi{x, si) X 
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ds2 G^(x, si + S2) exp(-S2/so) 
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(15) 
(16) 
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A similar computation can be performed for the third 
cumulant, 



III. ELEMENTARY PHENOMENOLOGY OF 
"POISSON STRING" 



r)3 r + oo 

{Jf)^A8csG^-^^ilsl dsGKx,s). (17) 

^OS "'-CXJ 

One can further notice that for parity reasons (7!) = 0. 
For the first component of the shear we however have 

{Jl)^8c,soG'§^^.l^, (18) 
^os ^ 

{j!)=-48c,.slG^§^^f^. (19) 

One can notice that the dimcnsionlcss skewness 

JdL^^^fE^y^'ji. (20) 

^^2)3/2 VS^IJ ,3/2- ^'^) 

The ratio c^jSj"^ being a priori finite (and surely the 
one-point probability distribution function of [i has no 
reason to obey Gaussian statistics), it implies that the 
reduced skewness vanishes at distance much larger than 
the correlation length sg from the string. This is a nat- 
ural expectation since at finite distance and for sq — ^ 
such quantities as the shear are obtained as sums of in- 
finite number of sources. So, very generically, whatever 
may be the statistical properties of /x, we expect that lo- 
cally, all quantities follow Gaussian statistics as soon as 

In the following we will thus assume to be in this 
regime, that is we assume that the energy distribution 
along the string has a vanishing coherence length. In this 
case the exponential factor can be replaced by a Dirac 
delta function so that the 2-point correlation function of 
/i can be written, 

(/i(si)^(s2)) = 2 C2 So Mo '5Dirac(si - S2) (21) 

and the other correlation functions are then not relevant 
for the description of the distortion properties. (It does 
not mean that the probability distribution function of /i 
is assumed to be Gaussian distributed!). Obviously, this 
result does not depend on the details of the small scale 
correlation properties of /i. One can see that in this limit 
only one extra parameter is required to describe /i, that 
is C2So- Since it is not a priori possible to distinguish 
between the coherence length so and C2, in the following 
we use Eq. ( ^l]) assuming that C2 = 1/2. 

Together with Eq. (||) we call this model of strings, 
the Poisson string model. We are aware of the dramatic 
simplifications we have adopted to describe the outcome 
of an a priori very complicated physics. Still we think 
that some of the basic properties of lens phenomenology 
can be captured by such a model in a much more real- 
istic way than a simple straight string with a uniform 
distribution. 




FIG. 1. Numerical experiment showing the amplification 
map, i.e. 1/ det(^), of a "Poisson string". The brightest pix- 
els correspond to infinite magnification: they form the critical 
lines. The darkest pixels correspond to a magnification close 
to zero. The solid lines correspond to the caustics, positions 
of the critical lines in the source plane. The external dashed 
lines are the counter images of the critical lines. 

On Fig. |l] we depict an example of numerical imple- 
mentation of such a cosmic string, showing various fea- 
tures associated with this lens system. It is obviously 
beyond the scope of this article to present the general 
theory of gravitational lenses (see |l^). The aim of this 
section is to give some hints on the expected phenomena 
encountered for such a model of energy distribution. 

The grey levels show the variation of amplification 
given by the determinant of A^^. Along the brightest 
areas the amplification is infinite; these locations form 
the critical lines. This is where the most dramatic lens 
effects could be detected: giant arcs, merging of images... 
Note that such lines simply do not exist for strings with 
uniform density! The critical lines are mainly made of 
two long lines running along the string without crossing 
it. This latter property is due to the fact that whatever 
the local value of /i the displacement field has always an 
infinite gradient at the string position. This point will 
be quantified more precisely in the next section. 
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Such a system of critical lines is associated with caustic 
lines in the source plane. These caustic lines are shown 
on Fig. H as thick solid lines: the left one for instance is 
obtained by applying the displacement field to the right 
side long critical line. If a background object happens 
to lie on this line, it will appear as a highly deformed 
object (with formally infinite magnification) on the im- 
age plane. The caustic lines are also directly related to 
the possibility of having multiple image of background 
objects. Each time one is crossing the caustic lines the 
number of associated images is changed by 2: at large 
distance the number of image is one; after one crossing it 
has 3 (in other words if an object is located in between 
the two caustic lines, it has 3 images), after 2 crossings 
it has either 1 again or 5 (this is the case for the caustic 
lines associated to the inner critical lines), etc... 

Thus, in general, on the image plane the critical lines 
are betrayed by dramatic lens effects such as merging of 
two images of background objects. It is then worth noting 
that we are here generically in a regime of 3 images in the 
vicinity of the string (except in rare cases where it can 
be 5) instead of 2 as for a strictly uniform string. The 
dashed lines show the locations of the counter images 
of the critical lines. It delimits the region, in the image 
plane, within which multiple images can be found. It can 
be noted however that the amplification rapidly decreases 
in the vicinity of the string, so that central images (i.e. 
the ones situated in between the two infinite critical lines) 
are expected to be strongly de-amplified, except when 
they are close to one of the critical lines. 

Moreover, the examination of Fig. ^ shows that nu- 
merous small images generically appear along the string. 
The presence of these images are due to the fluctuating 
small scale structure of the string. They are associated 
with an infinite number of critical lines (and caustics) 
near the string, that are only partially exhibited on Fig. 
due to the finite resolution of our simulation. In a re- 
alistic string these images would appear down to scales 
corresponding to the coherence scale sq. One aim of the 
next section is again to quantify more precisely this as- 
pect. 

We have checked that, in the limit of the resolution 
of this simulation, the global feature we see here do not 
depend on the shape of the one-point probability distri- 
bution function of /i neither on the actual resolution. 

IV. STATISTICAL PROPERTIES 

The aim of this section is to quantify some aspects of 
the results described in the previous part. Indeed, in 
general systems of critical lines depend non-trivially on 
the lens strength and also on the optical bench config- 
uration. So, at this point, it is not clear how typical is 
the plot presented on Fig. |^, and whether it depends or 
not on the parameter choice. Other aspects that may be 
difficult to grasp with a simple numerical experiment are 



the mathematical properties of the critical lines. Since 
the strings has an infinite number of substructures (when 
the coherence length is assumed to vanish) the properties 
of the critical and caustic lines may not be numerically 
stable. In particular there seem to be a large number of 
small critical lines close to the string. So we do not know 
a priori whether those lines form a fractal like structure, 
or whether their total length is finite. 

In the following we want in particular to locate the po- 
sition of the critical lines, compute their average length 
(per unit string length) as well as the length of their 
counterparts in the source plane, the caustic lines. And 
quite surprisingly, what in general would be an impos- 
sible task (because critical lines form complex patterns 
that depend nonlinearly and nonlocally on the energy dis- 
tribution) turns out to be computable to a large extent 
for this model. 

The calculations are primarily based on the fact that 
all local quantities obey Gaussian statistics and can thus 
be entirely described by their second moments. This 
property is valid because we assume to be in a regime of 
an infinitely small coherence length on the string. This 
would not have been the case otherwise, and not only 
should the correlation effects of n along the string be 
taken into account but also the fact that ^(s) itself is not 
a priori Gaussian distributed. 

The calculations of the moments of the displacement, 
the shear components and the shear gradients are the 
required ingredients and are given first. 

A. Correlation matrices 

1. Displacement field 

Not surprisingly we can check that, 

(ei(x)) --eosign(a;) (22) 

(Ux)) = (23) 

where 

eo = 47rG^A^o, (24) 

which corresponds to the case of a uniform cosmic string. 
We introduce the angular variable that defines the 
average energy scale of the string. The second moment 
of the displacement field is given by, 

Stt 1 

{^^{x)^J{x)) = —So ^il 5rj = Co ^tj- (25) 

One can note that the displacement fluctuations have the 
same amplitude in directions along the string and per- 
pendicular to the string (we have no convincing physical 
interpretation for that). A consequence of these results 
is that angular pair separations is no more fixed, (and 
given by 2^0); they may fluctuate along the string, and 
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the amplitude of the fluctuations is all the more large 
that the images are close to the string. The image pairs 
may also not be strictly orthogonal to the string direc- 
tion. The amplitude of these fluctuations depend on the 
dimensionless ratio sq/^o- The larger it is, the larger the 
relative fluctuations are. 

Such effects can be observed on Figs. |[ where pairs 
can be clearly identified, but distances and orientations 
indeed fluctuate from pair to pair. 



2. Amplification matrix 

As we already stressed, the main novel effect compared 
to a uniform straight string is the emergence of complex 
distortion features. The latter are related to the value 
of the local shear. For a uniform string the shear com- 
ponents are always zero. A reminder of this property is 
given by the average values of the elements of the ampli- 
fication matrix that all vanish. 

The shear components are explicitly given by 



71 



D 

72 = -16G 



OS 

D 



iy 



LS 



OS 



ds fi{s 



{y - s)x 



{{y-sr + x^y 



(26) 
(27) 



Their second moments can be calculated^. They depend 
on the distance to the string and are given by 



{7zix)-/jix)) = 



(28) 



1 



(32) 



and the critical lines are the locations where 7=1. From 
the relation (|30|), we can compute the probability that 
a point at a given location is within the critical zone 
(7 > 1) or out (7 < 1). This probability is given by. 



'Pcrit.(a;) 



7'(7)d7 = exp 



1 



exp 



2 so 4' 



(33) 



This probability has a very simple shape. It reaches unity 
nearby the string. This expresses the fact that the string 
itself is always in the critical region: no critical lines 
can actually cross the string, and there must always ex- 
ist two infinite critical lines running on each side of the 
string. The situation observed on Fig. |l] is thus com- 
pletely generic in this respect. 

This expression also indicates the behavior of typical 
critical line distance to the string. It must scale some- 
how as (so^o)^^"^ since this is the only distance that in- 
tervenes in Eq. (|33|). The exact calculation of such a 
quantity cannot however be obtained from the mere dis- 
tribution of the shear. The number of intersection points 
between the critical lines and any horizontal lines for in- 
stance is also not contained in the probability. Somehow 
the spatial correlations of the shear have to be taken into 
account. We will see in the following that the statistical 
properties of the shear gradients allow the computation 
of such quantities. 



For completeness we can also compute the cross- 
correlation terms, 



3. Shear gradients 



1 



Att X' 



■ So ?0 ^^3 ' 



(29) 



with the displacement field. These results allow us to 
gain some insights into this system. For instance one can 
compute that the local shear distribution function. It is 
given by, 



7'(7)d7 = exp 



7 



7d7 



where 



ill) -ill)- 



(30) 



(31) 



Since the local convergence vanishes, the local ampliflca- 
tion ^J, — 1/ det A is given by 



^It is actually possible to also compute analytically the two- 
point correlation functions, (^'yi{x,y)'yj{x' ,y')'^, but it has no 
use in the following. 



They are a priori four different shear gradient com- 
ponents. However, since the local convergence is zero, 
there are simple relations that relate these components 
together. 



dxl2 = dyjl, 

dyj2 = -dxii, 



(34) 
(35) 



so that only two quantities have to be considered. We can 
furthermore define a pseudo- vector 5 as 5 = {dx^i, dyji). 
The calculations of its statistical properties can be done 
in a similar way as previously, and one gets 



{Si{x)Sj{x)) = — so^l (5y, 
3 

{-f,{x)5j{x)) = ^^:^^0^0^Z3, 



(36) 
(37) 



from which we deduce that the correlation coefficient r 
with the shear field is given by 



5 



{12S2) 



^/3 



-0.86. 



(38) 



As a result we can explicitly write the joint probability 
distribution of the local shear and shear gradients, 



exp 



2(1 -r2) I ct2 



5^ 2rj-d 



-I 



(39) 



where r is given by (38), is defined in Eq. ( |3l| ) and 
as is given by = = {5l). 

B. Properties of the critical lines 

1. Where are the critical lines? 




0.5 



0.0 
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FIG. 2. Shape of the function ricrit.(a;) as a function of the 
distance from the string in units of (sq^q)^'^^- 

As mentionned previously we are now in position to 
address issues related to the critical and caustic lines. 
First thing that can be considered is the number den- 
sity, rTcrit.l^;), of intersection point between the critical 
line and any horizontal axis. This function will allow us 
to compute not only the typical number of such inter- 
section points but also the average distance between the 
critical lines and the string, and its fluctuations. The 
method we introduce here will also be employed in the 
next paragraph to derive the total length of the critical 
lines. 

If Xi denote the locations where det(^) = (which is 
equivalent to 7 = 1 when k = 0) on a given horizontal 
axis, the local number density of intersecting points is 



^crit 



{x) = ^ foiiac(a; - Xi), 



(40) 



where each Xi is an implicit function of the random vari- 
ables /x(s). Such a quantity is actually a random quantity 
that depends in particular on the y position. But what 
we really want is the average number density of these 
points, that would be measured as a geometrical average 
over y position, 



dy n„it.{x), 



(41) 



in the limit of a large L. If we assume this system to 
be ergodic, then the geometrical average of Eq. ( |4l| ) can 
be replaced by an ensemble average over the positions Xi 
(that are random quantities depending on the variables, 
7, (5, etc..) as, 

'^crit. {x) — 

To perform such an average, one needs to make a change 
of variable from Xi to local parameters 7 and 5. This 
is possible since we need to do this change of variable 
only in the vicinity of Xi (we follow here the same line of 
computation as what has been done in 3D or 2D Gaus- 
sian fields for computations of number density of maxima 
[Q). The constraint {x — Xi) — Q can then be replaced 
by a constraint on the random variables 7 and 5. Since 
at first order in [x — x;). 



7(0;) I + [x - Xi) dxl{Xi) 

we obtain that 

«crit.(a;) = 

7'(7,5V7d'^^^Dirac(7-l)|9.7l, 



(43) 



with 



(44) 



(45) 



The integral (^4|) can then be easily computed explicitly 
using the expression (p3). It leads to 



n„;t.[x) = V^e" 



(46) 



- TTx Erf(W — a;2)e 



when X is taken in units of (sq Co)^^"^- shown on Fig. || 
where we have depicted rlcrit. {x), the intersections mostly 
take place at a distance of about 0.75 (so Co)^^'^- The av- 
erage number of intersection points (on one side of string) 
is then 



H-oo 



da;ncrit.(a;) = 



1.155. 



(47) 
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The number of intersection points being an odd number, 
it means that the critical Hne crosses one horizontal line 
more than once in at most 7% of the cases. It supports 
the fact that we are dominated by the two long critical 
lines located on each side of the string. In rare cases, 
inner critical lines can give rise to complex multiple image 
systems. 

The typical distance of the critical lines to the string 
can also be computed. It is given by 



/o^°°dxa;ncrit.(a;) 



/q dxncrit.(a;) 
whereas the scatter of this distance is about 



0.70 (so Co')'/' (48) 



(49) 



One recovers the scaling behavior suggested in previous 
paragraph. Remarkably the average distance and the 
scatter of the distance follow the same scaling law. It 
suggests that the critical line system is somehow univer- 
sal, in the sense that it does not depend on the dimen- 
sionless ratio sq/^o- This idea is further supported in the 
next paragraph where we compute the total length of the 
critical lines. 



2. Length of the critical lines 

The length of the critical lines (per unit string length 
i) is a priori given by 

/>oc />L 

icrit. = 2 / dx dySBiiac{x-Xc,y-yc) (50) 
Jo Jo 

where (xc, yc) describes the running position of the crit- 
ical lines in the (x, y) plane and the factor 2 accounts to 
the fact that we take into account both side of the string. 
To complete this calculation one can perform a change 
of variables from (x, y) to (uc,Vc) where Uc and Vc are 
the coordinates of a point in a basis given by (n, t) , nor- 
mal and tangential vectors to the critical line at {xc,yc) 
position. Then the integral (M) reads 



dUcdVc (5Dirac(Uc)- 



(51) 



As previously the (^Dirac function can be replaced by a 
constraint on 7. In this case the modulus of the gradient 
along the x direction is replaced by the one along the n 
direction which is | V7I by construction. Finally one gets, 
returning to the x and y variables, 



Lcrit. = 2 / dx dy V{^,6)d'^d'6 x 
Jo Jo J 

|V7|foirac(7" !)• (52) 

This can also be computed straightforwardly by noting 
that 



so that 



IV7I = 75. 



72(51) 



(53) 



(54) 



The integral over y is then straightforward and gives a 
factor L. The integral over S can be performed without 
too much difficulty and gives. 



/ dxx^/2/o 
Jo 



exp 



57rx'^ 



(55) 



where Jq is the Bessel function. Finally, using an integral 
representation of Iq, one finds 



V3 



E 



L w 2.80 L, 



(56) 



where E is the complete elliptic integral. Remarkably 
this coefficient is independent of the parameter of the 
model. In particular it does not depend on the dimen- 
sionless ratio sq/^q. Because this integral is finite, it also 
proves that the closed critical lines have a finite total 
length despite the fact that they are in infinite number. 
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FIG. 3. From the same numerical experiment we show the 
variation of the location of the critical (right side) and caustic 
(left side) lines when the source plane is displaced. The grey 
levels show here the amplitude of the shear field. We see 
that the critical lines tend to stick around points where 7 
vanishes (and the points where it happens are independent on 
the source plane distance). The thin lines show the variation 
of the position of the counter image of the caustic lines. 

C. Properties of the caustic lines 

If one wants to detect such a string with distortion 
effects, the efficiency of such a detection will be directly 
proportional to the length of the caustic lines which are 
the locations on the source plane where the background 
objects are significantly distorted. 

The calculation of the length of the caustic lines is a 
priori a much more difficult task since it is related to 
quantities in the source plane. One can actually circum- 
vent this problem. The idea is that a length clement 
d^crit. on the critical line is going to be transformed in 
a length element on the caustic line dZcaus. according to 
the amplification matrix linear transformation 

d^caus. = A ■ d/crit.- (57) 

We are interested here only in the change of length |d/crit. | 
to |d^caus.|- It depends locally on the eigenvalues and 
eigen-directions of A. On the critical line the eigenvalues 
of A are and 2 (one is necessarily by construction, 
the other is given by the trace of the matrix, that is 2, 
since k = everywhere). As a result we deduce that 

Idlaus.l =2|dlcrit. -621 (58) 

where dZcrit. • 62 is the component of dZcrit. on the eigen- 
direction associated with the eigenvalue 2. The eigenvec- 
tor associated with this eigenvalue is 

- ( ) ■ 

If we set 

7 = 7 (cos -i/", sin-;/;) (60) 

and 

V7 = 7(5 (cos (/), sin(/)) (61) 

we deduce that 

^ \cos{cj)-^)+cos^\. (62) 

Id/crit.l v/2(l+7i) 

The shear field has no prefered direction and thus it is 
possible to integrate over ip with a flat distribution to get 



I d^caus. I 

|d/crit.| 

/ o |cos(0- V) +COS-0I = - (63) 

J-^ 27r ^(1 + 71) TT 

(this can be obtained by the change of variable, ip — 2t). 
Remarkably this result is independent of (p so that we 
directly conclude that 

4 

-^caus. = - icrit. ~ 3.56 L. (64) 
TT 

^From an observational point of view we have obtained 
with this result some insights into the detectability of 
such a string from large distortion effects. The cross- 
section of a Poisson string of length L can be estimated to 
be about 3.6 L times typical galaxy sizes, if galaxies are 
the objects that one uses to reveal the distortion effects. 

D. Cusps number density 

One could pursue these investigations by computing 
the number density of cusps along the critical lines, that 
is points corresponding to the merging of 3 images. These 
points are characterized by the property that the criti- 
cal line is locally along the eigenvector associated with 
the eigenvalue equal to zero. Technically it corresponds 
to the cases where |d/caus. | vanishes in previous calcula- 
tions. Since the detailed calculation of such a quantity 
requires one further derivatives of the shear, we do not 
present it here. We can however argue from simple scal- 
ing arguments that 

ncusps oc (soCo)"^^^- (65) 



E. Scaling laws 

Note that for our description to be valid, the typical 
distance of the critical line to the string, (sq^o)^^'^! must 
be larger than the coherence length along the string sq 
(assuming C2 to be finite). This is the case when the 
average displacement is larger than sq. 

Changing the position of the source plane is equiva- 
lent to changing the amplitude of while keeping sq 
constant. It acts indeed as a normalization factor (see 
Eq. |2^). If this amplitude is multiplied by A, the angular 
distance of the critical line to the string scales like, 

rfcrit. OC (66) 

But still the total length of the critical lines remains un- 
changed. On Fig. 3 we illustrate this effect by changing 
the amplitude of ^o- We see that the critical lines (thick 
lines on the right side) are drifting away from the string. 
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However they tend to get stuck on points where the lo- 
cal shear vanishes (dark patches on the figure) creating 
invaginations towards the string position. When gets 
large enough the critical line is disconnected and leaves a 
small inner closed critical line around 7 = points. This 
effect can be observed in particular in the central part of 
the figure. As a result the critical line has larger position 
fluctuations, but over a larger scale along the string, so 
that its total length is kept constant. One consequence of 
this mechanism is that there are locations on the image 
plane (close to points where 7 = 0) which correspond to 
quasi-superposition of critical lines from many different 
source planes. That should make them more likely to be 
detected. 

At the same time the caustic lines are sweeping over 
the source plane; their distance to the critical line, and 
to the string, increases like A, if the typical displacement 
is larger than (sq^o)^^'^- In this case we can also observe 
that the size of the multiple image region is determined 
by the displacement amplitude and thus scales like A (po- 
sitions of the thin lines on the right side). 

V. DISCUSSIONS 

This investigation provides a new description of the 
phenomenological properties of lens physics for cosmic 
strings. The model of "Poisson string" we present is 
based on general results which state that the lens effects 
of string are those obtained from lineic energy density. 
We think such a model must have captured most of the 
generic properties expected for string lens phenomenol- 
ogy as the resemblance with the results of numerical ex- 
periments (see Fig. 3 of |l^) strongly suggests it. We 
are however aware that the validity of the description we 
adopted has not been demonstrated (numerical simula- 
tions of string networks anyway lack the necessary reso- 
lution for describing correctly the small scale structures 
of strings). It is clear for instance that the projected po- 
sition of a string on the sky does not follow a straight 
line. We have checked that the results we found are not 
strongly affected if a non-zero curvature radius is intro- 
duced, as long as this radius is larger than the other 
quantities involved in this problem, i.e. ^0 and (so^o)^^'^- 

The results we have obtained suggest novel strategies 
for detecting strings: multiple images are still present, 
but images can have undergone significant distortions; 
small images of background objects can be found in large 
numbers... The most adapted strategy depends however 
on the resolution of the imaging apparatus and on the 
energy scale of the string. On Fig. ^ we show an example 
of an HST deformed image. In this case the energy scale 
of the string expressed in units of ^0 is similar to the 
angular scale of a galaxy. For GUT strings, <^o would be of 
the order of a few arcsecs which is indeed a typical galaxy 
size in deep surveys. The probability of observing such 
an event depends obviously on the survey coverage. The 



expected angular length of strings to be present within 
redshift unity has been estimated in and is such that 
on average one string is expected to cross a 100 deg^ 
survey. But even though such surveys are within reach 
in the near future, the probability with which one can 
detect a string (or conversely the rejection one can put 
on the existence of strings with a given energy scale) still 
deserves further investigations in particular on the effects 
of multi source planes (the amplitude of the lens effects 
depends on the position of the source plane) that may 
partly blur the searched effects. 
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FIG. 4. Example of a deformed image by a Poisson string. The field corresponds to an external region of cluster A2218 
(taken by the Hubble Space Telescope). If such a field was put at 2 = 1, then an intercepting string at z ~ 0.8 along the 
line-of-sight would produce multiple images as observed on the 2 pictures (the typical pair separation is about 5 arcsec). The 
right panel shows this system with a resolution of about 0.5 arcsec, whereas the left panel shows the same image with sub 0.1 
arcsec resolution. In this case the distortions of the background galaxies are clearly exhibited. Note also that a number of 
small images appear along the string. In high resolution images, that might be the most effective way of detecting strings. 
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